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1. Introduction 

Many years ago [1-4] we suggested a picture of the formation and composition of 
the QCD string as a linear chain of gluons, which are the perturbative excitations 
the theory. This picture is motivated by ’t Hooft’s large-A c expansion. In partic¬ 
ular, a time-slice of a high-order planar diagram for a Wilson loop (Fig. [I]) reveals 
a sequence of gluons, each of which interacts only with its nearest neighbors in the 
diagram by mainly attractive forces. This immediately leads to the idea that the 
QCD string is composed of a “chain” of constituent gluons, each held in place by its 
attraction to its two nearest neighbors in the chain. The challenge in such a model 
is to understand how the attractive force between gluons, which is essentially due 
to one-gluon exchange, can manage to hold these massless constituent particles in 
a bound state. In this article we address this problem, and in particular we show 
how recent perturbative results for the static quark potential, and the related force 
renormalization scheme [5], bear on this issue. We propose an appealing way to 
interpret perturbative QCD that leads to a self-consistent extrapolation of the per¬ 
turbative force of a quark on an anti-quark to a linear confining force. We do not 
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Figure 1: The gluon chain as a time slice of a planar diagram (shown here in double-line 
notation). A solid hemisphere indicates a quark color index, open hemisphere an antiquark 
color index. 

claim that this in any way proves confinement, but rather that it provides a model 
framework for thinking about the physics of confinement which stays conceptually 
close to perturbation theory. A computational scheme is presented, involving vari¬ 
ational and perturbative elements, which we hope to eventually apply to calculate 
the ground state (and the corresponding static quark potential) of the gluon chain. 
Even without a full calculation of this ground state, we are able to use our variational 
framework to demonstrate both the logarithmic broadening of the QCD flux tube, 
and the existence of a Liischer —c/R term at long distances. 

Section 2, below, reviews the motivation of the gluon-chain model, with emphasis 
on how this model accounts for features of the confining force which are problematic 
for other theories of quark confinement. In section 3 we discuss in more detail the 
physics of gluon chain formation, and the binding of gluons in the chain. The running 
coupling is, of course, crucial to this dynamics, and at the coupling strengths relevant 
to chain formation there are important issues of renormalization scheme dependence 
that must be confronted. In section 4 we present a variational framework for semi- 
perturbative calculation, and show how it applies to a string tension calculation. In 
this section we demonstrate how the gluon-chain model accounts for roughening, and 
the Liischer — c/R term in the static quark potential. The last section contains some 
concluding remarks. 

2. The Gluon-Chain Model 

Like any theory of quark confinement, the gluon-chain model aims at explaining the 
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linearity of the static quark potential. However, it is now widely recognized that 
in addition to the linearity feature, there are at least three other properties of the 
confining potential which a satisfactory theory of confinement is obligated to explain: 

• Casimir Scaling: Consider the potential between static quarks in a repre¬ 
sentation r of the gauge group. From the onset of linearity in the potential, to 
a finite (adjoint string-breaking) scale, the string tension a r is proportional to 
the quadratic Casimir C r of the representation, i.e. 

o> = a F (2.1) 

L-f 

where the subscript F denotes the fundamental representation. For rectangular 
LxT Wilson loops with L/T fixed, the range of L for which the Casimir scaling 
law is valid increases logarithmically with N c . 

• Center Dependence: Asymptotically, the string tension can depend only on 
the N-ality of the group representation r, i.e. on its transformation properties 
under the center subgroup of the gauge group. 

• String Behavior: The diameter of the color-electric flux tube between static 
sources is believed to grow logarithmically with the separation L of the sources 
(roughening), and there is a — c/L contribution to the asymptotic potential 
(the Luscher term) which is due to quantum fluctuations of the QCD string, 
rather than Coulomb attraction of the quarks. 

Taken together, this is a challenging set of conditions. The abelian monopole 
theory [6], for example, has a very hard time accounting for Casimir scaling [7], 
as well as for the center dependence of certain operators [8]. Instanton [9] and 
meron [10] mechanisms are consistent with Casimir scaling, but have difficulties with 
center dependence. The center vortex theory [11] is in perfect accord with center 
dependence, and it is at least roughly compatible with Casimir scaling, as shown 
in ref. [12]. On the other hand, the vortex theory does not really explain the high 
degree of accuracy of the Casimir scaling rule, which has been revealed in numerical 
simulations [13]. Finally, the string-like behavior of the QCD flux tube seems to pose 
problems for any theory of confinement based essentially on one-gluon exchange (e.g. 
the proposal of ref. [14]), as well as for the proposal of stochastic confinement [15]. 

The gluon-chain model of QCD string formation, on the other hand, meets the 
conditions listed above in a rather simple and appealing way, as we will discuss below. 

2.1 Linear Potential 

As a heavy quark-antiquark pair move apart, and their color charge separation in¬ 
creases, we expect that at some point the interaction energy increases rapidly due to 
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Figure 2: As heavy quarks (large circles) separate, energy is minimized by keeping the 
average color charge separation below a certain limit. This is achieved by pulling out 
a sequence of gluons (small circles) between the heavy quarks. Again, solid and open 
shadings denote quark and antiquark color indices, while dotted lines indicate contracted 
color indices between a quark and a gluon or between neighboring gluons. 

the running coupling (cf. section 3). Eventually it becomes energetically favorable 
to reduce the effective charge separation by inserting a gluon between the quarks. In 
the N c —> oo limit, the quark and antiquark can only interact with the intermediate 
gluon, but not directly with each other. As the heavy quarks continue to move apart, 
the process repeats, and we end up with a chain of gluons, as shown in Fig. |2j. in 
which the average distance between color charges remains fixed, irrespective of the 
separation of the heavy sources. The energy of the system is approximately NE g i uon , 
where N is the number of gluons in the chain, and E g i uon is the kinetic and nearest- 
neighbor interaction energy per gluon. If the quarks are separated by a distance L, 
and the number of gluons per unit quark separation ( N/L = 1 /R) is fixed, then 

E chain « NE gluon = = aL (2.2) 

where a = E g i uon /R is the string tension. The linear growth in the number of 
constituent gluons is the origin of the linear potential in the gluon-chain model. 

Alternatively, we may understand the linear potential in terms of a constant 
force. The force between neighboring constituent gluons is dependent on the average 
separation R = L/N of gluons along the quark-antiquark axis. If this separation 
remains fixed as L increases, then the average intergluon force remains fixed. The 
intergluon force, which is the same everywhere in the chain, can be interpreted as a 
string tension, which is constant irrespective of the quark separation. 

2.2 Casimir Scaling 

To leading order in N c , a group character in representation r is given by a product 
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Figure 3: Adjoint string-breaking in the gluon chain model. Two gluons in separate 
chains (I) scatter by a contact interaction, resulting in the re-arrangement of color indices 
indicated in II. This corresponds to chains starting and ending on the same heavy source. 
The chains then contract down to smaller “gluelumps” (III). 

of group characters in the fundamental representation 

Xr[g\ oc (xf[. 9]) + sub-leading terms. (2.3) 

By factorization at large- N c , a Wilson loop in representation r has a string tension 

cr r = M r ap (2-4) 

at N c —> oo, where M r = n + n. In this limit, the quadratic Casimir is C r = M r N c j 2. 
Exact Casimir scaling is therefore a property of the planar limit. 

The gluon-chain model, which is motivated by large-iV c considerations, inherits 
this property. A heavy source in representation r is the terminus of M r separate gluon 
chains, one for each of the n quark and n antiquark charges in the direct product 
forming the representation r. Since the chains do not interact in the N c — oo limit 
(the interaction is a non-planar process, as can easily be verified by considering the 
relevant Feynman diagrams), the total energy is simply the sum of the energies of 
each of the chains. In this way, Casimir scaling is obtained, at least at large N c } 

2.3 Center Dependence 

Asymptotically, the string tension of static sources in a higher color group repre¬ 
sentation r depends only on the N-ality, rather than the quadratic Casimir, of the 

-‘-At small N c , the Casimir is not simply proportional to M r . On the other hand, at small N c , 
interactions between the chains cannot be neglected. We cannot say, at present, whether this model 
predicts any substantial deviation from Casimir scaling at small N c . 
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group representation. Let us consider how this comes about in, e.g., the adjoint rep¬ 
resentation (the analysis is easily generalized to other representations). Beginning 
with heavy quark sources in the adjoint representation, there are two gluon chains, 
as shown in Fig. [|. Nearby gluons in each chain can scatter, e.g. by a contact in¬ 
teraction, and rearrange the sequence of colors as shown in the figure. The result 
is that two gluon chains transform to two “gluelumps”; the adjoint string is bro¬ 
ken (providing the sum of gluelump masses is less than that of the double chain), 
and the resulting string tension is zero. This is the correct prediction, since the 
N-ality of the adjoint representation is also zero. However, the scattering process 
indicated is non-planar, and the transition rate from the two-chain structure to the 
gluelump structure is 1 /iV 2 suppressed. In the large- N c limit, we therefore recover 
exact factorization and Casimir scaling. 

2.4 String Behavior 

It should be obvious, just from the figures, that a gluon-chain is a discretized string 2 
of some kind, with the constituent gluons playing the role of “string-bits.” Therefore 
it is reasonable to expect that, due to quantum fluctuations of the chain configura¬ 
tion, we should find the logarithmic broadening of the color-electric flux tube with 
quark separation (roughening), as well as the Liischer —c/R term term in the static 
quark potential. These effects are very non-trivial, however, and need to be demon¬ 
strated in the context of the gluon-chain model. We will postpone the analysis to 
section 4. 


3. The Force Renormalization Scheme and Gluon Chain For¬ 
mation 

In the remarks above, we have passed lightly over a fundamental issue. Gluons are 
massless particles. The question is how nearest neighbor forces between constituent 
gluons actually manage to bind such particles together in a chain. Clearly it is 
hopeless to find a binding mechanism at very weak coupling a s = g 2 /47r, since the 
kinetic energy of a gluon confined to a region of size R is of order 1 /R, whereas the 
interaction energy (at tree level) is only of order —a s /R. The fact that the effective 
a s is not really constant, but grows with distance scale i?, is obviously of crucial 
importance. 

Superficially, a coupling which grows monotonically with distance seems exactly 
what is needed for quark confinement. However, different choices for what the cou¬ 
pling measures lead to drastically different estimates of the quark-antiquark static 
potential. For example, if we say that a s (R ) is a measure of the potential energy 

2 For a treatment of discretized bosonic strings in light-cone gauge, see ref. [16]. 
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between quark and anti-quark (the “V-scheme”), we would then write the familiar 
expression for the potential energy of a static quark and antiquark separated by 
distance R in an SU(iV c ) gauge theory 




V V?/ 2 R 


(3.1) 


It should be understood that this formula is a definition of the running coupling, and 
there is a complicated relation (beyond two loops) of off. ( R ) in the V-scheme to the 
running coupling in, e.g., the MS scheme [5,17]. 

If al ( R ) grows monotonically with R, as suggested by the first few terms in the 
renormalization group Gell-Mann-Low function, the potential in eq. ( |3.1| ) is actually 
the opposite of what is required for confinement. Assuming that af(R) grows with 
i?, perhaps even blowing up at some finite R^ (the Landau singularity), the potential 
in eq. ( |3.1|) leads to a force which first weakens, then vanishes at some point, and 
finally becomes ever more strongly repulsive. The problem is that a monotonically 
increasing af (/?,) tends to drive the potential away from the V = 0 axis in the 
direction of negative V, whereas in fact the static potential crosses the axis at some 
point and becomes positive. Therefore, as defined in the V-scheme, af(R) derived 
from the true static potential cannot possibly grow monotonically. Instead, at some 
scale, af must start to become smaller and eventually change sign. There seems to 
be little hope of relating such behavior to the running coupling, at least up to three 
loops in perturbation theory. 

As long advocated by one of us (C.B.T.), it is less problematic to define the 
running coupling by the force (the derivative of the potential) that the static quark 
exerts on the static anti-quark [18] 



N c a s (R) 
2 R 2 


( _ J_\ vrACR) 

V K) 2 i? 2 


(3.2) 


where we have defined the ’t Hooft coupling held fixed in the N c —> oo limit by 
A = N c cx s /7r. In the rest of this paper a s and A will always refer to the running 
coupling defined through the static force. This scheme, also advocated in recent 
years by Sommer [19], has the obvious advantage that the running coupling so defined 
doesn’t have to change sign. The corresponding Gell-Mann-Low function, defined as 

m = (3.3) 

is given to three loops for Nj = 0 but for any N c by [5] 

A) = -La 2 - La 3 - (3.795...)A“ - ■ ■ • (3.4) 

The exact coefficient in the last term is rather cumbersome, so we have merely quoted 
its numerical value to 3 decimal places. 
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Figure 4: The force between static quarks in units of the Sommer scale Rq ~ 0.5 fm, 
computed perturbatively to three loops at N c = 3 in the force renormalization scheme [5]. 
Also shown is the asymptotic string tension of (430MeV) 2 x Rq. 


Calculating A (R) by truncating 0 at three loops leads to a force that, as R 
increases, drops to a minimum at R m given by 


2A (Rm) + 0(A (R m )) — 0, 


(3.5) 


after which the force increases , blowing up as (R^ — R)~^ 3 at the Landau singularity 
R = Roo, as shown (for N c = 3) in Fig. [|. In this figure, Rq is the Sommer scale 
Ro ~ 0.5 fm. The corresponding static quark potential, in the force renormalization 
scheme, is obtained by integrating the force, 


V q ,(R) 


*R 


Vgq(R A ) 


dR 


1 Ra 


F(R) 


(3.6) 


where Ra -C Rq and V^Ra) may be estimated at sufficiently small a s (RA ) by using 
eq. ( |3.1|) , with a]i « a s . Note that the Landau singularity in the force is integrable 
as R Roo from below, and the integrand becomes complex for R > R^. Thus the 
potential curve stops at R^ at a finite value. It was found by Necco and Sommer in 
ref. [5] that the resulting three-loop potential is a surprisingly accurate match to the 
lattice Monte Carlo result, almost up to the Landau point at Roo/Ro ~ 0.78. This 
three-loop potential (with R A = 0.15f? 0 ) is displayed in Fig. [3]. 
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Figure 5: The static quark potential at N c = 3, computed perturbatively to three loops 
in the force renormalization scheme. 


It is clear that a potential of this sort can trap a massless particle, whose kinetic 
energy |p| decreases as 1/R by the uncertainty principle: Minimizing 1/R + V(R) 
gives 


1 

& 


V'(R) 


(,_i_ \zm 

\ m) 2 & 


(3.7) 


in other words A (R) = 9/Att ps 0.71 for N c = 3. The effective coupling required for 
binding a massless particle is of order 0(1), but not huge. There is thus hope that, 
after taking binding into account with such an extrapolation of the RG improved 
perturbative force, subsequent corrections may be under control. It follows that, 
according to the force renormalization scheme, it is possible to have bound states 
containing massless constituent gluons, such as glueballs and gluon chains. Moreover, 
there is no indication of any gross failure of perturbation theory up to R ~ R m . On 
the contrary, as just mentioned, it seems that perturbation theory fits the numerical 
lattice data quite well. 


3.1 Saturation Mechanism 

Beyond R = R, rn the three loop force increases, and the corresponding potential is 
concave upwards. In contrast, a confining force is expected to be asymptotically 
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constant: One can say that the perturbative force extrapolated to R > R m over¬ 
confines. In fact, there is a rigorous theorem, analogous to the well known convexity 
of the thermal free energy, which states that the static force never increases with 
quark separation, or, equivalently, that the static potential is concave downwards [20]. 
It appears that physics beyond three loops, at R > R m , must supply a mechanism 
to weaken the extrapolated perturbative force law, not strengthen it. In our opinion, 
this mechanism is supplied by the formation of a gluon chain, as described in section 
2. As the quark-antiquark separation increases, it can become energetically favorable 
for a gluon to materialize between the quark and anti-quark. Because the gluon is 
in the adjoint representation, its color can be arranged so that it is simultaneously 
a sink of the color flux from the quark and from the anti-quark, and indeed this 
is precisely the arrangement of color dictated by the N c —> oo limit. The effect 
of a gluon with this color orientation is to shield the direct force of the quark on 
the antiquark, and thereby reduce the effective separation of color charge. Further 
increase of R leads to the successive appearance of more constituent gluons, creating 
a gluon chain between the quark and anti-quark. Thus the separation of directly 
interacting color charges never exceeds some maximum value. It is presumably this 
saturation mechanism which prevents the static force from eventually increasing with 
quark-antiquark separation, in accordance with the concavity theorem. The same 
saturation mechanism which bounds the average color charge separation also bounds 
the corresponding effective coupling. If this maximal effective coupling is not too 
large, it should be possible to treat the interaction between neighboring gluons by 
perturbative methods. 

In the next section we develop a variational approach to the calculation of the 
properties of the gluon chain responsible for quenching the increasing force beyond 
R m . But it is worth noting here that the form of the perturbative force law shown in 
Fig. 0, plus the hypothesis of a saturation mechanism, already provides the basis for 
a zero parameter estimate of the string tension, namely a ~ F(R m ). In other words, 
we accept the extrapolated perturbative force up to the point where it violates the 
concavity theorem, after which we replace it with the simplest behavior consistent 
with the theorem, namely a constant. This reasoning is reminiscent of the Maxwell 
constuction in statistical physics, which restores convexity to an approximate calcu¬ 
lation of the free energy of a Van der Waals gas. 

As already noted, there is a theorem which tells us that the increasing part of 
the force curve in Fig. [| (at R > R m ) cannot apply to the actual force between static 
quarks. However, if that curve is taken as a paradigm of the nearest-neighbor force 
binding gluons in the chain (neglecting change in gluon number), then its increasing 
part may play a role in providing a restoring force towards an equilibrium gluon 
separation, particularly if the mean intergluon separation in the chain exceeds R m . 
In this case, there is a net restoring force on any gluon fluctuating away from the 
mean separation. Again, the saturation mechanism will prevent any overall increase 
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of the static force with quark separation. 

To arrive at an estimate of the string tension, it remains to find A (R m ) and R m 
from Eq. |3T5| . Truncating 0 at three loops leads to a cubic equation for X(R m ), with 
the numerical solution A (Rm) ~ 0.540, so (1 — l/N^)nX(R m )/2 ~ 0.754 for N c = 3. 
This value of A gives the relative size of the terms entering into Eq. |3.5| as 

2A - 0 O A 2 - 0iA 3 - 02A 4 = .540(2 - 0.990 - 0.413 - 0.598), (3.8) 

so it is clear that one is extrapolating beyond the strict validity of perturbation 
theory. However, if the three loop term turns out to be anomalously large, the 
prospects for our perturbative approach would be brighter. From the slope of 
the meson Regge trajectories a' ~ 0.86GeV -2 we have the estimate from data of 
a = l/2na' ~ (430MeV) 2 . This leads to R m ~ (495MeV) _1 ~ 0.40fm. According to 
Figures 2 and 4 in ref. [5], we can infer that with Aj^g = 238MeV, the phenomenolog¬ 
ically preferred value, R m ~ 0.3fm. The string tension for this value would be higher 
by a factor of (1.33) 2 ~ 1.76, i.e. a ~ (570MeV) 2 . Quark production, vetoed in our 
N c —> oo model, is expected to reduce the predicted tension somewhat, but probably 
not this much. Conversely, the string tension inferred from the meson trajectories 
would require Aj^g = 177 MeV. 

Since the saturation mechanism must come into play at R ~ R mi where pertur¬ 
bation theory begins to fail, we would conclude that the QCD ground state in 

the presence of static qq sources is dominated, at R ~ R m , by a one-constituent gluon 
component (perturbation theory expands around a state with zero constituent glu¬ 
ons). 3 If that is so, then the average separation of color charge along the axis joining 
the static sources will be R m /2. Assuming this “axial separation” of color charge re¬ 
mains fixed as quark separation L increases, as it should according to the gluon chain 
model, then for large separations, at N c = 3, we may estimate L/N ps R m /2 sa 0.15 
fm, where N is the number of gluons in the chain. This does not mean that the 
gluon is sharply localized on the axis, however. Indeed, in order to be bound it must 
be spread to a size r at which the running coupling is of 0(1), which means r > R m . 


4. A Variational Framework, and String-like Behavior 


In order to make the notion of “constituent gluon” precise, it is necessary to work in 
a Hamiltonian framework, which also means picking a physical gauge, e.g. Coulomb 
or light-cone gauge. The ground-state wavefunctional in any physical gauge can be 
expressed in the path-integral form 


T 0 [A(x)] 


DA^t < 0) 5 F{A ) 


A^p[A] exp 



(4.1) 


3 An old lattice Monte Carlo investigation of the gluon content of the QCD string (second article 
of ref. [2]) suggests that the one-constituent gluon component dominates beyond 3 lattice spacings 

at (3 = 2.4, corresponding to R = 0.73 R 0 . 
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where La is the gauge-held Lagrangian, F(A) = 0 the gauge condition, and A pp the 
Faddeev-Popov determinant. 4 An excited state can be constructed by multiplying 
the ground state with some polynomial in the helds, i.e. 


ty ex [A\ = Q[A}^ 0 [A}. 


(4.2) 


For example, the wavefunctional appropriate for a glueball state would be 

Nmax 

= E («) 

N=1 

where is the component of the glueball state with N constituent gluons, 

which has the general form 



dxidx 2 ...rfxjv / WM 2 ... AtJV (xi,x 2 , ...,xjv) 
TrA w (x 1 )A /12 (x 2 )...i4 /1JV (x JV ) 1 v h 0 [A], 


(4.4) 


For our purposes, a “constituent gluon” simply refers to a gluon held operator mul¬ 
tiplying the exact ground state. The energy of any excited state ( |4.2| ) above the 
ground-state energy, is given by 


(ggxggex) 

{^ex\^ex) 


(H). 


(4.5) 


Defining 


Qt = Q[A{x,t)] (4.6) 

the energy of the excited state can be computed from a correlation function in the 
Euclidean-time version of the theory 


E 


\Hlr log(Q ' TQ - T) 


(4.7) 


where ( Q' t Q_t) is the Euclidean vacuum expectation value of these operators. In 
the perturbative calculation of E one need only include the connected diagrams. 

Equation (4.7) is the starting point for a held-theoretic variational approach to 
bound states of constituent quanta [3,21], The idea is to hrst make some trial ansatz 
for Q[A] involving a few parameters, calculate the energy E tria i of the trial state 
perturbatively, via eq. (O, and then find the parameter values which minimize 
this trial energy. The result is a variational estimate for bound-state energy, and an 
approximation to the corresponding wavefunctional. 

4 We note that in a physical gauge either App is trivial (as in light-cone gauge), so that there is 
no need to introduce ghost fields, or else the ghost fields do not propagate in time, as in Coulomb 





In the case of the gluon chain, we take the static sources to be at points x = 0 
and x = L, and consider ^-constituent gluon states generated from trial operators 
of the form 


Qckain[A] = </ ai (°)| I dxidx 2 ...dx X VWa...|«(Xl, x 2 ,...,Xiv) 

^r i ( x i)' 4 r s (x2)...^r' ;+, (x«)y"'(L) <4.s) 


where ^ is a trial N-gluon wavefunction. In the N c —> oo limit, only interactions 
between nearest-neighbor gluons (in the diagrams) need be taken into account. A 
consistent chain picture will require that the RG-improved intergluon interaction 
energy as a function of gluon separation is qualitatively similar to the static quark 
potential shown in Fig. Whether this happens remains to be seen. 

4.1 A simplified model 

We will not attempt, in this article, a variational computation of the gluon-chain 
energy £ in the full field theory. The perturbative evaluation of (Q^Q-t) for the 
gluon chain state is very tedious, especially beyond tree level, and again involves 
delicate issues of renormalization scheme dependence. However, we can capture 
most of the important qualitative features of the field theory calculation with a 
simple quantum-mechanical model, in which the (N-l)-gluon Hamiltonian is taken 
to be 

N -1 

H = J^lpn 


N -1 


n =1 


) r(x„ - x„_i) + Us(xi) + V qg { L - X/V_i) 


(4,9) 


n =2 


where V(;x) is the gluon-gluon interaction energy, while V qg (x) is the interaction 
energy between a static quark and its neighboring gluon. The restriction to nearest- 
neighbor interactions is, of course, justified in the large N c limit. In field theory 
the interaction between each pair of neighbors is described by renormalization group 
improved Bethe-Salpeter style ladder exchanges. In addition to the RG improved 
Coulomb exchanges that build up the static force displayed in Fig. [| exchange of 
transverse gluons and contact interactions will add spin and momentum dependence 
to the interaction. We can anticipate that the most attractive spin channels are those 
in which neighbor gluons are in a spin 0 state, so that for a chain we should have 
anti-ferromagnetic order [22], As we would do in the field theory calculation, we use 
the variational method in this simplified model to obtain an approximation to the 
N-gluon ground state energy £(L). As already mentioned, if the potentials V and 
V qg have the qualitative behavior shown in Fig. ||, a self-consistent gluon chain will 
arise. 
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• Product Ansatz — Relative Coordinates 


We begin with a simple product ansatz for the gluon-chain wavefunctiona.l 


< P(xi,x 2 ,...,x A r_i) = A JjV'fUi 


(4.10) 


where the {uj} are relative coordinates 


Uj = X,; - Xj_i 


(4.11) 


x 0 = 0 , x^v = L 


(4.12) 


It is convenient to also change integration variables to relative coordinates, and the 
integration measure becomes 


d 3 q N ] / 

dV = ——- d 3 Ui exp iq ■ ( ^ u, - L 

' ' i= 1 


dV |^| 2 = 1 (4.13) 


where integration over q gives a delta, function enforcing the constraint 


Ui = L = JVR 


(4,14) 


Defining 


F( q) = / d 3 u\'ip(u)\ 2 e 


2 jq u 


(4.15) 


and the Fourier transform of the relative-coordinate wavefunction 

r d^i i 


(4.16) 


we find for the normalization constant 


V 2 = / T^e-'i ^fq)" 


(4.17) 


and the expectation value of the gluon nearest-neighbor potential 


2 f [ d 3 u\^{ri)\ 2 V (u)e iq ' u 


(V(uk)) = A ‘^ F(q) 


(4.18) 
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For the purpose of computing the expectation value of kinetic energy, it is useful 
to introduce 


(e iyPk ) = A 2 


= A 2 


(<KP k - p)> = A 2 


d3q e-^F(q) N - 2 


( 27 r) 


d 3 ud 3 u'e iq - (u+u '^* (u)^*(u')^(u + y)i/>(u' - y) 
d3q e~ iqL F(q) N ~ 2 


(2vr 

^ d 3 kd 3 k / e iy ' (k ~ k ' ) (/)*(k + q)0*(k' + q)0(k)0(k') (4.19) 

d3 d 


r e-* qL F(q ) JV - 2 


(2vr) 3 

[ d 3 k(j)*( k + q - p)</>*(k + q)0(k - p)0(k). 


(4.20) 


In the limit IV —» oo we can evaluate the g integral by a saddle-point technique. Call 
q 0 the saddle point determined by 

^(qo)=ih = i R . (4.21) 

Then we have in this limit 

(V(Ufc)) = F(qo)- 1 J d 3 u\^(n)\ 2 V(n)e iqoU (4.22) 

f(u)f(u^(u + Y)^(u' - Y) (4.23) 

(|Pfc|) = F(q 0 )~ 2 J d 3 p\p\ I ^e- jY p I d 3 u J d 3 u’e iq °< u+ 
ip*( u)ip*(u')ip(u + Y)ip(u r — Y) 

= F(q 0 )~ 2 J d 3 p\p\ J d 3 k 

<t >*(k + q 0 - p)0*(k + q o )0(k - p)0(k). (4.24) 


Note that in our product ansatz for the trial wavefunction, all dependence on the 
gluon number in the chain, indicated by the subscript k, has disappeared in V(iifc) 
and |pfc|. The total energy of the trial state is then simply 5 

s = W«|p*|> + (C(u fc ))). (4.25) 

5 We disregard the distinction between V and V qg at the ends of the chain, on the grounds that 
this is unimportant for £/L at large L. 
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Once this quantity is minimized as a function of the variational parameters in the 
wavefunction, our estimate for the QCD string tension is the energy per unit length 

= £ = ilM±h>h) . (4,26) 

L R 

To go further, we need to choose a definite ?/;(u), f° r example the gaussian 

^(u) = e -“ 2/2r2 . (4.27) 


In this case there are two variational parameters. One of them is the parameter r in 
the above gaussian. The other is the number N of gluons in a chain between heavy 
sources separated by distance L or, equivalently, the distance R = L/N. For the 
gaussian wavepaeket we find 


0(k) = r 3 e~ k2r2/2 
F( q) = (nr 2 ) 3/2 e- q2r2/ \ 

— t 

Rr^/n Jo 


qo = 


<V) = 


2iR 

r 2 

^2 /„2 


uduV{u) ( ! r - e ~ iu+R) /r 


for central O(u). In particular we find 
|u|/ RryJH 


du (VO-«) 2 A 2 - e - (u+R)2/r2 ) = 


erf (R/r) 

R 


(4.28) 


(4.29) 


To evaluate the kinetic energy, note that 


(k + q 0 - p) 2 + (k + q 0 ) 2 + (k - p) 2 + (k) 2 
= 4k 2 + 4k • (q 0 - p) + (q 0 - p) 2 + q^ + p 2 
= (2k + q 0 - p) 2 + qo + P 2 


so the k integral is a simple Gaussian leading to 


2 \ 3/2 


(IPfcl) = ( I dp\p\e 


-rV/2 _ i / £ 


r V 7r 


Thus the product ansatz with gaussian wavepackets leads to 

£ 1 / 8,1 /tA 

L ~ rR \ n + R^' 

In particular, taking for V the instantaneous Coulomb potential 

Qtq 


V <“> = - C ^l 


we have 


£ 


1 18 CfOLi 


rR\7r 


R 2 


erf 


(R\ 
\r J 


(4.30) 

(4.31) 

(4.32) 

(4.33) 

(4.34) 
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• String Wavefunction Ansatz 


The second type of trial state we consider here is the ground state wavefunction 
of a discretized string. Strings with discrete degrees of freedom were studied some 
time ago in ref. [16], and we will borrow some results directly from that reference. 

The ground state of the discrete string is the state annihilated by all the lowering 
operators for string modes 


«; )( vl ; ( x i.x 2 ,....x.v- i) - 0 j A 1 ( 4 -35) 

where the position and momentum of the fc-th gluon are related to the string modes 
via 


x * = N k 


Pfc = 


2 1 



NTr 


N -1 I - 

( a - - a U sin icY k ) 

m= 1 V 


E 




f a, n + a m 


) (w* 


where 


com = 2 sin ■ 


rmr 

2N 


is the frequency of the mth mode. Then we have 

(0|e'*' Pfe |0) = exp(-x 2 Aj) 

1 \ 3/2 

exp 


where 


(0|<5(p — Pfc)|0) = 


pM 


4 A \) 


(0|e ip ' Ufc |0) = exp (iR ■ p - p 2 A*) 

3/2 

k I eX P 


(0|5(u - u fc )|0) = 


1 


47tA? 


u — R) 

4AJ 


N-l 


A\' = 


NTn 


m= 1 


E . TTL7T 9 I ,nn 1 1 

sm —— cos —— \ k H— 
2 N \ N \ 9 


rri'K 

TV 


4TVT 0 
1 


7r 1 7 t(4/c + 3) 1 7t(4/c + 1) 

cot-1— cot-cot- 

AN 2 AN 2 4iV 


An = 


7tT 0 

0 


N 


E 


(4/c + 3)(Ak + 1) J ttTq 
T 0 


rp N-l 

in \ . rmr . 2 /rmr, 

sm-sm - k 

2N V N 


m= 1 


7T 


. Za 

4A 

To 

_ 

7T 


7T 1 7t(4/c + 1) 1 7t(4A; — 1) 

cot-cot-1— cot- 

AN 2 AN 2 AN 


1 + 


(4k + l)(Ak — 1) J 7r 


T 


o 


(4.36) 

(4.37) 


(4.38) 

(4.39) 


(4.40) 

(4.41) 

(4.42) 

(4.43) 
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From these expressions we find, for V{u) proportional to l/|it| as in eq. ( |4.33|) , 


£ 

L 


1 4 
R 0T 


V^2 


Cpa s 

R 2 


erf 



(4.44) 


Having tried two different types of trial wavefnnctions, the question is which 
leads to a lower ground state energy. The answer is that the string wavefunction is 
the better of the two, at least for inter-gluon potentials proportional to l/|u|. To see 
this, simply define r = 2y/Ai, so that 


£ 

L 


1 8 
rR s/tt 


a/ A1A2 


CpOls 

R 2 


erf 



(4.45) 


and, for the string wavefunction, A 1 A 2 = 1/n 2 . Comparison with eq. ( |4.34| ) shows 
that for any values of r and i?, the string wavefunction has a slightly lower energy 
than the relative-coordinate product wavefunction. 


4.2 String Tension 


Without making any assumptions about the inter-gluon potential, the energy per 
unit length of the gluon-chain in the string-wavefunction ansatz is 


£ 

L 


1 r 8 1 

R 7T 3 / 2 r 


<V»> 


(4.46) 


As explained in section 2, the minimal value of this energy per unit quark separa¬ 
tion, as a function of the parameters r, R , is the variational estimate of the QCD 
string tension a in our quantum-mechanical model. Of course, even the simplified 
model requires as input the intergluon potential H(u), and to get this interaction 
energy right one should probably use the variational approach in the full field theory. 
However, even without knowing the intergluon potential precisely, we can make some 
rough estimates using knowledge of the three-loop static Coulomb potential as input. 
First of all, if V{u) were simply due to tree-level effects, then we would have 

(V(u))oc/A\. (4.47) 

Likewise, the running coupling is generally regarded as a function of the inverse color 
charge separation [5,17]. Let us therefore assume that (V) depends on the variational 
parameters r, R only through the 

1 

s 

so that we may write 


VEV 
1 


= -2-erf (- 
u| / R \ r 


V(s) = (V(u)) 
. dV 
F (S) S dF 


(4.48) 


(4.49) 
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where F(s) can be regarded as the magnitude of a semi-classical “force” between 
gluons. 

Minimizing S/L results in two conditions. The first, obtained from minimizing 
wrt R , is 


or, in view of (|4.46j) , 


I | H- 1 

+ 

( 5 ) 

) - M™ 

(4.50) 

- 

a™ 

(4.51) 


£ 

L 

This equation has a simple physical interpretation. On the one hand, the QCD string 
tension in the gluon chain model is simply the energy of the chain per unit quark- 
antiquark separation, i.e. £/L on the lhs of (|4.51|) . On the other hand, the “tension” 
in this system should be related to the change in potential energy of the system with 
respect to small deformations 5R in inter-gluon separation; this is the “restoring 
force” in the system along the quark-antiquark axis. The minimum condition ( |4.51| ) 
equates these two types of expression for the string tension. The second condition 
for the minimum is obtained from minimizing £/L wrt r, which gives 


1 = 


7J-3/2 j,2 Qj, 


(4.52) 


Solving (|4.52| ) for F(s) and inserting into (|4.51| ) gives an expression for the string 
tension (er = £ / L) 


ds /<9s\ 


-i 


where 


dR \dr ) 7r 3 / 2 r 2 


= erf l (—)-^=—erf 2 (R/r) exp(— R 2 /r 2 ) 

oR r \/tt r 


(4.53) 


ds 

dr 


2 R 2 


7T r 


erf (R/r) exp(— R 2 /r 2 ) 


(4.54) 


Given F(s), we could use (|4.51|) and (|4.52|) to determine R and r, and then ( |4.53|) 
would give the string tension. Of course, this procedure assumes that eqs. ( |4.51| ) and 
( |4.52|) have a solution at finite r, R. For F(s) determined at tree level, without taking 
account of loop corrections or the running coupling constant, eq. ( |4.52| ) at fixed R 
will typically have a solution at r = oo (small a s ), or r oc R (large a s ). But then 
( |4-51|) will only be solved at R = oo or R = 0 respectively. Clearly, incorporating 
the running coupling is a crucial ingredient to a self-consistent chain solution. Even 
though the detailed strengths and shapes will be different, we can expect the RG 
improved forces (with -0 truncated at one, two or three loops), corresponding to 
V and V qg in attractive channels, to have the same qualitative behavior as Fig. fj]: 
an initial decrease to a nonzero minimum, and a blowup at finite R due to the 
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Landau singularity. This qualitative behavior is sufficient to guarantee a consistent 
chain solution. Quantitative accuracy relies on the chance that the effective coupling 
required for the solution is not too large. 

At tree-level, the inter-gluon force F(s) has three components: the Coulomb 
force, the magnetic exchange force, and the contact interaction. Some previous 
variational [3] and bag [22] calculations suggest that the magnetic exchange force 
between neighboring gluons with total spin 0 will be attractive, and roughly equal in 
magnitude to the Coulomb force, while the contact interaction is repulsive and (for 
glueball wavefunctions) about 60% as large as the Coulomb force. For an accurate 
estimate of all contributions, we must actually carry out the calculation for the 
gluon-chain by the field theory method outlined above. However, based on the work 
in refs. [3,22], we expect that the sum of contributions results in a (renormalization- 
group improved) intergluon force of the form 


F(s) w (4.55) 

s 1 

where k is a number which is probably between one and two. 6 

Even without knowing F(s) precisely, we can still use the three-loop perturbative 
results for the static potential to make some educated guesses about the size of 
R and r, and arrive at a “ballpark estimate” for a. At the end of section 3, we 
argued for R = R m /2 m 0.3 Rq , on the grounds that a constituent gluon should 
appear between the quark-antiquark pair when they are L = R m apart, since this is 
where the perturbative expansion around a zero constituent gluon state seems to be 
breaking down. For the r parameter, we just note that it is not reasonable to have s 
beyond the Landau singularity, because before that distance we would expect more 
constituent gluons to have appeared and reduced the effective charge separation. On 
the other hand, the bound state condition will require a running coupling a s (s) which 
is 0(1), otherwise the gluon kinetic energy dominates, and there is no binding. So s 
cannot be much less than the Landau singularity. A guess which is perhaps not so 
far off is s — R m = 0.6i?o which, when combined with R = R m ,/2 = 0.3i?o, implies 
r = 0.63i? 0 . These particular guesses happen to land close to the right answer, i.e. 

1.26 
a -J% 

where the accepted answer for the QCD string tension is a = (430 Mev) 2 = 1.18/Aq. 
This rather close agreement is probably fortuitous, since we are only making guesses 
for r and R, but it does show that our numbers are in the right ballpark, and there 
is some hope that the eventual calculation of cr will be in reasonable agreement with 
the phenomenological value. It is also worth noting that with these parameters, the 

6 At loop level, there will also be corrections to the kinetic energy term, which tend to increase 
the size of that term [3]. 


(4.56) 
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spread in the wavepacket of each constituent gluon (set by s), is roughly twice the 
size of the average intergluon separation along the quark-antiquark axis which is 
R 0.3i?o- Thus the QCD string is rather thick even before roughening effects come 
into play. 

Continuing a little further, we note that the above choice of parameters R , r 
implies 


a(s) 


J_i e R 2 y 
kCf vt 
1.60 
kC p 


(4.57) 


If we use the 3-loop running coupling in ref. [5], and N c = 3, then s ~ 0.6i? 0 would 
require k ~ 2, which is probably a little large. But none of these numbers should 
be taken very seriously at this stage, and a full field-theoretic calculation using the 
string wavefunction ansatz is clearly required, before we can draw any quantitative 
conclusions. 


4.3 Roughening and the Liischer Term 

The phenomenon of roughening, i. e. the logarithmic growth with L of the transverse 
size of the gluon chain, requires the long range correlations contained in the string 
wave function ansatz. For example, the product ansatz would not display this effect. 
Indeed it is easy to calculate 


<0|xL|0> = 


rs-/ 


D- 2 
2NT 0 

D — 2 
2ttT 0 


^-4 sm 2 (rmrk/N) 
sm(mn/2N) 


m 


In N — 


r 2 (D — 2) 



(4.58) 

(4.59) 


These same correlations also produce certain finite size effects, that are sub-dominant 
for L/R —> oo, in the trial energy. These come from explicit 1/iV corrections and 
from the fact that 2 have mild k dependence. 

To estimate these effects we first note some sums: 


N -1 




To 

2 


N -1 

Emu 

sin —— 

9 N 


m= 1 




To 

2 


2 N 

71 


1 

2 


71 

24N + "' 




fc=l 


1 

2 To 


N-l 

mn 
2 N 

m= 1 



rs_/ 


1 

2 n 


2 N 

7r 


1 

2 


71 

24iV + "' 


(4.60) 

(4.61) 


It is of interest to estimate the finite size contributions to the energy with our string 
trial wave function. For this purpose the detailed form of the kinetic and potential 
energy is largely irrelevant and it is more efficient to use an arbitrary kinetic K( p) 
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and potential V (r) energy in onr formulas. We then have, after changing integration 
variables to dimensionless ones: 


(K(pk)) = 


(V(u k )) = 


3/2 


7 r 


- I d 3 pe~ p2 K (2pyfAl 
I d 3 ue~ u2 V (r + 2uy / Af 


(4.62) 

(4.63) 


Then calling the limiting values of A* 2 —> Ai ;2 , we can expand the energy to first 
order in A( : 2 — Ai i2 , and do the sums over k. 


N-l 


N 


£ = E< a '<p + 


k =1 


k =1 


= (JV -1) 


7T 3 / 2 

y-y A 2 — A 2 \J a 2 


I d 3 p e- p2 K (2 Pv /A^) + N -ij J d 3 u e~ u2 V (r + 2x1^ 


fc=i 

v 

E 

fc=i 


A 2 7T 3 / 2 
Af — Ai v^Ai 

A^ "^T 2 " 


d 3 p e p2 p • VA” 


P 


d 3 u e- u “u • W 


(N - 1)-^ J d 3 p e~ p2 K (2py^j + N 
2 1 VA~2 r 


71 


3/2 


(4.64) 

d 3 u e~ u ~V ("r + 2u-\/Ai N ) 


+ 

+ 


7T 7T 

~~ 4 ~ 48iV 


7r 7T 

’4 ~~ 48iV 


7T 3 / 2 

2 n I 


7r 


3/2 


I d 3 p e P “p • V/l 
d 3 « e" u2 u • W. 


In these expressions V is always the gradient with respect to the argument of the 
function that follows it. Next we minimize S/L at fixed L with respect to the 
parameters T 0 and R = L/N, and obtain 


0 = 


V^2 

7T 3 / 2 


d 3 p e' pl p • V/b - 


v^Ai 

7T 3 / 2 


d 3 u e _u u • VI/ 


(4.65) 


° T 2 


J d 3 p e p2 K ^2py/A^ + J d 3 u e U "V ^R + 2ua/A^ 

+ 1 j d 3 u e - u2 ^-W (R + 2u v^aT’ 


(4.66) 


These two equations determine To and R. When both are satisfied, we can simplify 
our expression for the variational energy: 


£ = NR- 


3/2 


7r 


J d 3 u e" u2 ^ • W (r + 2 Uv /A^ --L jd 3 p e~ p2 K {ipsfK^ 
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+ 1 “ I" 2 In ^J d3p e ” p2p ' va '( 2 p\/^) ( 4 - 67 ) 

= ^J d3u e ’ u2 f’ W ( R + 2uv/ ^) _ |S^/^ e ' p2 P - VX(2p\/A^) 

+ I 1 - \\ / d3p e " p2p • VAr - / % e ^ K ( 2p v^) • ( 4 - 68 ) 

The first two terms are the linear potential and Liischer terms respectively. The 
last two terms are independent of L and can be absorbed in the source energy. For 
A '(p) = IpI, as appropriate for massless gluons, the integrals in the last three terms 
can be immediately carried out, yielding 

£ = ^f d3u e ” u2 f' w ( r + WaT) 

7t AR ' 7 r 4 

24L rypK - 2. r7r 3 / 2 

For the bosonic string the factor multiplying n/2AL in the Luscher term would be 
D — 2 = 2 in 4 dimensional space time. With our estimate above of R ~ r/ 2 this 
coefficient is estimated to be 2/« 1.13. But these numbers are far too preliminary 
to arouse disappointment in such a discrepancy, especially since it is not even settled 
that the bosonic string result is correct for QCD. 

5. Conclusions 

We have developed our picture of a gluon chain into a viable calculational framework 
for the physics of quark confinement which stays close to perturbative ideas. The 
crucial observation is that the RG improved force law shows behavior which can be 
interpreted as over-confinement. This shows that perturbative physics might contain 
the germ of quark confinement which, in combination with the gluon chain saturation 
mechanism, leads to a concrete proposal for detailed calculations. 

We sketched a field theoretic variational approach to the gluon chain wave func¬ 
tion. For a preliminary estimate of the the string tension and other features of the 
QCD string, we studied a simplified model, which replaced field theoretic exchange 
interactions with a quantum mechanical potential. This exercise shows how the 
chain model incorporates the important effects of roughening and the Luscher term 
in the quark potential. In addition, the numerical estimate for the string tension was 
certainly in the right ballpark. 

The next stage of this project is to redo the variational treatment completely in 
the context of field theory. One issue left unresolved is whether the RG improvement 
of the interaction between neighbor gluons on the chain shows the same qualitative 
behavior as that of the RG improved static force. If it does, the consistency of our 


(4.69) 


23 



physical picture will be confirmed. Whether the detailed quantitative predictions of 
onr model will then turn out as accurate as we hope — say to 15 or 20 per cent — is 
not yet clear. 

We should finally note that the gluon chain model, which is formulated in terms 
of particle-like excitations, is not opposed to a description of quark confinement in 
terms of some special class of field configurations, e.g. center vortices or monopoles, 
which dominate the QCD vacuum. Instead, these approaches should be seen as com¬ 
plementing one another. The center dependence of the asymptotic string tension 
provides a good example. In terms of particle (gluon) excitations, we can easily un¬ 
derstand that this dependence is due to string-breaking by particle production, which 
results in color screening of the higher-representation heavy sources by constituent 
gluons. In the gluon chain model in particular, the process is illustrated in Fig. |3j. On 
the other hand, if one tries to explain the center dependence of large Wilson loops in 
terms of field fluctuations affecting the Wilson loop holonomy, then it is clear that 
the area law must be due exclusively to fluctuations in the loop holonomy among the 
center elements of the gauge group. This leads (perhaps inevitably) to a picture of 
the QCD vacuum as being dominated at large scales by center vortex configurations. 
Thus we have both a particle (chain breaking) and field (center vortex) explanation 
for the same phemonemon, namely, the N-ality dependence of the asymptotic string 
tension. These particle/field descriptions need not contradict one another; they are 
more likely to be dual descriptions of the same underlying physics. 

The attractive feature of the gluon chain model is that it offers a simple and con¬ 
cise account of so many features of the QCD confining potential: linearity, Casimir 
scaling (at large N c ), center dependence, roughening, and the Liischer term. Be¬ 
yond that, the model provides a promising framework for quantitative calculation of 
the string tension and, perhaps, low-lying masses. Whether these calculations are 
practical, and if so how the results compare with phenomenology, remains to be seen. 
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